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Abstract 

We propose an approach which enables one to obtain simultaneously the 
glueball mass and the gluon mass in the gauge-invariant way to shed new light 
on the mass gap problem in Yang-Mills theory. First, we point out that the 
Faddeev (Skyrme-Faddeev-Niemi) model can be induced through the gauge- 
invariant vacuum condensate of mass dimension two from SU(2) Yang-Mills 
theory. Second, we obtain the glueball mass spectrum by performing the col- 
lective coordinate quantization of the topological knot soliton in the Faddeev 
model. Third, we demonstrate that a relationship between the glueball mass 
and the gluon mass is obtained, since the gauge-invariant gluon mass is also 
induced from the relevant vacuum condensate. Finally, we determine physical 
values of two parameters in the Faddeev model and give an estimate of the 
relevant vacuum condensation in Yang-Mills theory. Our results indicate that 
the Faddeev model can play the role of a low-energy effective theory of the 
quantum SU(2) Yang-Mills theory. 

Key words: glueball mass, vacuum condensate, Faddeev model, knot soliton, mag- 
netic monopole, quark confinement, 

PACS: 12.38.Aw, 12.38.Lg 



^ E-mail 
^ E-mail 
^ El-mail 
^ E-mail 
^ E-mail 



kondokSf acuity. chiba-u.jp 
a . onoQgraduate . chiba-u . j p 
akihiro . shibataOkek . jp 
sinoharaSgraduate . chiba-u . jp 
tomOcuphd.nd. chiba-u.ac.jp 



1 Introduction 

It is widely accepted that the quantized Yang-Mills theory, i.e., quantum gluody- 
naniics (QGD) is the fundamental theory for the strong force. In the high-energy 
region, the perturbative expansion in the coupling constant gives an efficient method 
for calculations by virtue of the ultraviolet asymptotic freedom. In the low-energy 
region, however, some non-perturbative methods are needed to tame the strong cou- 
pling problem. It is rather difficult to perform the analytical calculation to derive 
non-perturbative results keeping the gauge invariance manifest. In such a case, it is 
sometimes useful to investigate the corresponding low-energy effective theory (LEET) 
instead of tackling the fundamental theory itself. 

In this paper, we demonstrate that the Faddeev (or Faddeev-Niemi) model PP 
can be used as a realistic LEET of SU(2) Yang-Mills theory by deriving the gauge- 
invariant glueball mass and gluon mass of the original Yang-Mills theory from the 
Faddeev model. In fact, it was suggested in a previous paper j2] that the Faddeev 
model is obtained as a low-energy effective theory of QGD by performing the non- 
linear change of variables (NLCV) of the gluon field ^ jSj and assuming the existence 
of novel gauge-invariant vacuum condensation (X^)ym of mass dimension two ^ |3]. 
The relationship between Yang-Mills theory and the Faddeev model is similar in spirit 
to that between QGD and the Skyrme model which expresses baryons as topological 
solitons of the meson field [Hj . 

The relevant NLGV for SU(2) gluon field ^(x) is represented by introducing the 
unit vector field n{x) = {ua^x)), (A = 1, 2, 3) as ^ 

£^^{x) =Cf,{x)n{x) +ig~'^[n{x),d^n{x)] +X^{x), 

2tT{n{x)n{x)) = 1, tr(n(a;)X^(x)) = 0, (1) 

where we have used the su{2) valued field: n(x) := n^(a;)T^, X^(x) := X^{x)Ta, 
Ta = |o"A; {<^A- Pauli matrices) and c^ and X^ are specified by n and is^ as 

c^{x) = 2tr(n(xX(x)), X^(x) = -ig~^[n{x), D^[s^]n{x)]. (2) 

The existence of gauge-invariant vacuum condensate of mass dimension two (X^)ym 
has the following implications. 

(1) The gauge-invariant gluon mass is deduced from the four-gluon interaction 
— ;^5f^(X^ X Xjy)^. For example, a simple self-consistent treatment in the mean field or 
Hartree-Fock appoximation leads to the replacement: 



{gX^ X X,) ■ ((7X^ X X'' 
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(3) 



^This vacuum condensation corresponds to a gauge-invariant version of the on-shell BRST- 
invariant vacuum condensation of mass dimension two in the modified Maximal Abchan gauge 
and the generahzed Lorentz gauge proposed in [S]. The existence of (X^)ym 7^ in the original 
Yang-Mills theory has been examined in numerical simulations 0] . The analytical treatment of this 
condensate will be given elsewhere and is not discussed in this paper. 

^This form of NLCV was known as the CFN decomposition |H1IZIE1 proposed by Cho [H] and 
Faddeev-Niemi [J] . In [Hj , three of the authors have given a new viewpoint of the CFN decomposition 
on which this paper is based. 



Therefore, X-gluons (off-diagonal gluons) acquire the mass Mx oc Jg"^ \^p) ^^'^ 

eventually decouple in the low-energy region below this scale. After integrating out 
the X-gluon field, the resulting LEET is expected to be written in terms of the vector 
field n{x) alone. This LEET will agree with the Faddeev model. ^ 

(2) The (gauge- invariant) glueball mass can be derived from the Faddeev model 
with the Lagrangian density: 



\Kl{d,nf - 1( 



^F = T^^lidf^nf - -C,[n ■ {d^n x 9,n)]^ (4) 



where the bare (tree) values of the parameters are given by [2] 



Af = v^(-X2)ym, C, = l/g'. (5) 

Then the two parameters Af and C4 in the Faddeev model are gauge invariant in 
the original Yang-Mills theory under this identification and the physical quantities 
calculated in the Faddeev model are automatically guaranteed to be gauge-invariant 
in the original Yang-Mills theory. 

(3) Quark confinement is expected to follow from the non-vanishing string ten- 
sion a (as a coefficient of the linear quark- antiquark potential V{r) = ar) which is 
proportional to (X^)ym: 

0- = const.(-Xj)YM- (6) 

There exists an intimate relationship between the glue-knot and the magnetic monopole 
which is a basic ingredient in the dual superconductor picture for quark confinement. 
This issue will be discussed in a subsequent paper. 

Thus, our approach presents an efficient framework for understanding the mass 
gap in Yang-Mills theory and quark confinement in QCD. In this paper, as a step 
along this line, we derive glueball mass spectra by performing the collective coordinate 
quantization of knot soliton solution in the Faddeev model. Glueballs are identified 
with knot solitons and their excitations. The knot soliton is a topological soliton 
with non-vanishing Hopf index Qh whose existence is suggested by the non-trivial 
Homotopy group 7r3(S'^) = Z. We can determine the two parameters Af and C4 of 
the Faddeev model by using two glueball masses Mo++ and M2++ for 0^^ and 2~^~^ 
(obtained by numerical simulations on a lattice) as inputs to predict glueball masses 
other than 0"*"^ and 2"'"+. Finally, we estimate the gluon mass Mx and estimate the 
vacuum condensation (— X^)ym of the original Yang-Mills theory. 

2 Classical knot soliton solution 

First, we obtain the knot soliton as a static and finite-energy solution n{x) of the 
Faddeev model with the energy given by 

E=J^^ d?x I^A^ (9,n)2 + ho^[n . {d,n x dku)]^^ (j, ^^ = 1, 2, 3). (7) 



^There are other approaches for deriving the Faddeev model from Yang-Mills theory. See e.g. 

[Tniini 



The energy E^j^^cAni^)] as a functional of n{x) for arbitrary values of the two pa- 
rameters Af and C4 obeys the scaling relation: 



1/2, 



-ir^i/2; 



EA,,cAn{x)] = AFCi^'E.^Ani^FCrx)]. (8) 

Once the solution n* minimizing E'l 1 is known, therefore, the solution corresponding 
to Sap,c4 is given by 



UAp^cAx) =n^{AFCi ' x), (9) 

and the mass of the knot soliton for arbitrary parameters obeys the relation 

Map,c4 = ^Ap,C4 Kp,C4] = ApCl^'^i,! [n*] = AfCI^'M,. (10) 

In what follows, * denotes the quantity calculated at Ap = 1 and C4 = 1. 
We adopt a simplified version ^2] of the toroidal ansatz J3] '■ 



n := (ni, 722,713) 



cos{mt, — nif), sm(7724^ — mp)^ 



(7n,7z = ±l,±2,---), 



Piv) + 1 



Piv) 



:iii 



where m and n are non-zero integers, and the toroidal coordinate {rj, ^, (f) is given by 



X =a 



sinh rj cos ip 



y = a- 



sinh 7] sin </? 



sin,^ 



cosh T] — cos C, cosh r/ — cos .^ 

< ?7 < CX), -TT < ,^ < TT, < (/? < 27r. 



cosh T] — cos ^ ' 



;i2) 



It is easy to understand that the Hopf topological charge of the Hopfion (static 
soliton solution with non-zero Hopf topological charges) under this ansatz is given by 



Q 



H 



mn 



(m,n 



±1,±2, 



(13) 



since the Hopf charge is equivalent to the linking number of two circles (in S*^) obtained 
as preimages of two distinct points in the target space S"^. See Fig. ^ 




Figure 1: The Hopfion configuration with Qh = 2 (m = l,n = 2) is indicated by a set of 
arrows on the torus specified by a toroidal coordinates. 



We consider the discrete and continuous symmetries of the color field n{x) under 
the toroidal Ansatz (lllj) characterized by indices n and m related to the Hopf charge 
Qh = nm. 

Planar reflection with respect to 



(a) yz plane: x -^ —x <^=^ if ^^ n — if 

ni(-x, y, z)m,n =(-l)"ni(x, y, z)^ 
n2{-x, y, z)m,n =(-l)"n2(x, y, z). 
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(14b) 




(14c) 



(b) xz plane: y 
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(c) xy plane: z 



ni{x, -y, z)m,n =ni{x, y, z)m-n = ni{x, y, z)_,n,n, 

n2{x, -y, Z)m,n =n2{x, y, z)m-n = -n2{x, y, Z)_rn,n, 
l^six, —y, Z)m,n =■'^3(2^) Z/) Z)m,n, 

ni{x, y, -z)m,n =ni{x, y, z)m-n = ni{x, y, z)_rn,n, 
n2{x, y, -z)m,n = - n2{x, y, z)m-n = n2{x, y, z)-m,n, 
iT-aix, y, —z)m,n ='n'3{x, y, Z)m,n, 



where n^^x, y, z)m,n = ri^i^x, y, z)-rn,n = n^i^x, y, z)m-n = n^i^x, y, z)-rn-n- 
Composing these two at a time yields the symmetries associated with 
Spacial rotations by vr radians about 



1 x-axis 



n(x, -y, 



<b)x(c)) 



-(y2n{x,y,z).m,n<y2 
-a2n{x,y,z),n,n(y2 



-e^'^'^/^nix, y, z),n,ne~'^''''''^ 
-e'^''^/^n{x, y, 2)^,^6"*^''^/^ 



(15a) 
(15b) 
(15c) 



(16a) 
(16b) 
(16c) 



[n = even) 
(n = odd) 
(17) 



(ii) y-axis (=(a)x(c)) 



n{-x,y, 



-a2n[x,y,z)m,nO-2 
-ain{x,y,z),ri,ncri 



-e^'^'^^^'^n{x,y,z)m,ne ^'^'^2/2 (7^ = even) 

.gi7rai/2^(^^^ y, Z)m,ne~'''''^^'^ {u = odd) 

(18) 



(iii) z-axis (=(a)x(b)) 



n — x, 



-y,z) 



'tT^yx, y, Zjjji^n 
a3n{x,y,z)m,n(^3 



[n = even) 
,i7ra3/2^(^2,^ ?/, 2)m,„e-^'''"='/2 (^ = odd) 

(19) 



where we have used (TjakCTj = —Ckij 7^ k),o-k{j = k) (no sum over j) following 
from CTjCrfc = 6jk'i- + iejk£(T£, and also the relation e^^^'^^^'^ = 1 cos | ± iuj sin | = 
±iaj. Making the product of all three rotations (i),(ii) and (ii) recovers the original 
configuration, as a trivial cross check. 



Finally, we consider Parity transformation. Making the product of all three refec- 
tions (a),(b) and (c) yields 
Parity: 

n{x, y, z)_m,n = -e''"'^/'^n{x, y, z)m-ne~''"'^^^ {n = even) 

Now we proceed to obtain the solution of the Faddeev model. After performing 
angular integrations over (f and ^ under the toroidal ansatz, we find that the energy 
functional reduces to 



/■oo 

E(Af, C4; a) =Ala / dr] 2 
Jo 



TT^ 



(u'Y 



U^ 



+ Kn,„(77)(l-u2) 



+ 6*40"^ / dr] 2TT'^iu' Y cosh. 7] siwh-qVmnW)-, (21) 

Jo 



where n^ = u{r]) := j^fhii^ u' := ^, and Vm,niv) ■= ^ + ?^^/sinh^?7. We obtain 
the knot soliton solution {n{x)} by varying the energy functional E = A^aE^"^^ + 
C4a~^E^^^ with respect to u and a, where we impose the boundary condition f{r] = 
0) = (X), i.e., the north-pole n = (0,0,1) at z-axis and /(?7 = 00) = 0, i.e., the 
south-pole n = (0, 0, —1) at the circle C := {x"^ + y"^ = a?,z = 0}. 

The scaling relation corresponding to (jHI is easily derived in the toroidal ansatz: 
E = AfC4 [dE^'^^ + d~^E^^^] using a dimensionless parameter d defined by 

a = Ap^d^^d. (22) 



The variation of E with respect to a yields a = J E^'^'i / E^"^^ which determines the 

radius a = Ap^Cl^^d = K^^Cl'"^ ^fs^^^JW) of the circle C, i.e., size of the knot 
soliton. This relation is substituted back into the energy functional to obtain E = 
2AfCJ v£'(2)^(4) ^hich is further varied with respect to u. The energy has the 
upper bound E < AfC^^E'^^^ + E^^^). 

In Fig. El we plot the numerical solution of the variable n-s = u{t]) := p^ )li for 
a knot soliton satisfying boundary conditions u{0) = 1 and u{oo) = —1. The energy 
density e{ri) as a function of rj is concentrated on a finite width which represents the 
location (thickness) of the knot. The rest mass M of the knot soliton is equal to the 
energy of the static knot soliton. Our numerical result (Table 1 below) shows e.g., 
M^ = 465 for m = l,n = 2;Qh = inn = 2, which is consistent with the previous 
results inilinilll. 



3 Quantization of knot soliton 

Second, we perform the collective coordinate quantization of the knot soliton solu- 
tion for the Faddeev model. The collective coordinate quantization for the Faddeev 
model has already been investigated in pT)|ll5j. However, the author of [THj has used 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
x=eta/(eta+1) 



Figure 2: The knot soliton solution u at Ap = 1 and C4 = 1 for m = l,n = 2; Qh = fnn = 2. 
(Left vertical axis) 713 = it vs. x = rj/{l + rf), (Right vertical axis) energy density e vs. x 
where /^^ (ir/e(7?) = ^(2) + ^(4)_ 



a different background configuration from ours. The authors of ^Q have given more 
general and mathematical framework for the quantization of the knot soliton than 
the collective coordinate quantization in the Faddeev model. Moreover, they have 
discussed the collective coordinate quantization too. However, the necessary expres- 
sions and the details of calculations to reproduce their results were not written in 
the paper ^Hj, which disables us from utilizing their results. In this paper, we have 
worked out all the steps of collective coordinate quantization to give explicit expres- 
sions and numerical data, following the way similar to jTTj. Furthermore, we give an 
argument for the relationship between the Faddeev model and the original Yang-Mills 
theory, although the paper ^H] aimed at a physical application due to the fermionic 
quantization of knot solitons which is quite different from ours based on the bosonic 
quantization. 

We observe that for a generic configuration n{x), there exists a nine-parameter 
set of configurations, all degenerate in energy (zero modes), obtained from n{x) by 
some combination of translation (X), spacial rotation (5), and isospin rotation {A): 



n(x] 



n [x] 



An{R{B){x - X))A\ A,B e SU{2) 



(23) 



where we have used the Lie-algebra su{2) valued field: n := uaTa, Ta = \o'a, and 
the spacial rotation is represented by orthogonal matrices Rjk{B) = ^ti^crjBakB^). 
The minimal or semi-classical quantization called the collective coordinate quantiza- 
tion proceeds by promoting the parameters X, A, B to the time-dependent dynamical 
variables X (t) , A(t) , B (t) to remove the classical degeneracy of the static configura- 
tion. Then the dynamical ansatz adopted in collective coordinate quantization is 

n(f) -^ n'(f , t) = A{t)n{R{B{t)){x - X{t)))A\t), A, Be SU{2), (24) 

where the prime serves distinguishing the dynamical variable n'{x,t) from the static 
background field n{x). The translational degrees of freedom X{t) can be ignored 
if the knot soliton is quantized in its rest frame. By inserting this ansatz into the 
Faddeev Lagrangian density ^p and integrating over the three-dimensional space. 



the Lagrangian in the rest frame is determined to be^ 

Lf= f d^x^F[n', dn'] =T-M, 
T =-ajUjkak + ajWjkbk + -bjVjkbk, (25) 

where the variables aj , bj. are defined by 

ttj := —itT{ajA'A) isospin rotation, bj := iti^ajBB^) spacial rotation, (26) 
and the inertial tensors are given by 

Ujk = / d^x [Al{5jk - rijUk) + C^diTijdenkj 

d^x {-2A2tr([r„ n][Tk, n]) + 2C4tr([[T,-, n], d,n][[Tk, n], d,n])} , 
Wjk = — d^x I Ap(ri X iLkn)j + CidinjU ■ {iLkfi x den) > , 
Vjk = / d^x JAp(zLjrl) ■ (iLkfi) + C^{iLjn x dfi) ■ {iLkU x din)\ , (27) 

with the angular momentum operator L = {Lj),Lj := ej^xuPe. = —i^jkt^kdi- Here, 
Wjk and Vjk are obtained from Ujk by the replacement: Wjk = Ujk[[Tk,'n] -^ Lku] 
and Vjk = Wjk[[Tj, n] -^ Lju] = Ujk[[Tj, n] -^ LjU; [Tk, n] -^ Lkn]. 
The scaling relation for the inertial tensor reads, e.g., 

V,k[n{S)]A,,c, = A^'Cl%k[n{A^'Cl^'x)],,,. (28) 

The inertial tensors under the toroidal ansatz have the following properties refiecting 
the symmetries of the classical background configurations considered in the above. 

i) Ujk, Wjk and Vjk are diagonals and all off-diagonal components vanish except 
for 1^21 = nWi2 ^ (n = ±1). 

ii) a) f/ii = U22, Vu = V22, 
b) V33 = nWss = n'^Uss . 

iii) Wii=0 = W22 . 

The property i) follows from the fact that the off-diagonal components j ^ k vanish 
after the angular integration over (p by its periodicity. The property ii-a) for the 
perpendicular components is a direct consequence of the cylindrical symmetry, while 
ii-b) for the parallel ones comes from the identify which holds for the toroidal ansatz: 

d 
L-sn = -i-K-n = -n[T3, n\. (29) 

The property iii) follows also from explicit calculations. The angular integration over 
^ is necessary to show Wn = = W22 for all integers n and m, since the angular 

■^Without the restriction to the rest frame, the kinetic energy T has extra terms of the type 
GjYjkXk and bjZjkXk, in addition to ^XjMjkX^- 



integration over y? alone leads to Wu = = W22 {n 7^ ±1) and Wu = —nW22 7^ 
(n = ±1). This is because the integrand is an odd function of C, where ^ -^ —C, 
corresponds to the reflection with respect to the z = plane, i.e., x -^ x, y -^ y, 
z — > —z. 

After lengthy calculations, we obtain the explicit form: 



a I T/fe 



Vxi =^n + Ki 



\u 



r\2 



V^^ = / dr] A'pa''n'{ coth^ Vr^-^ + 
Jo I 1 — n^ 



m^(3coth^r/-2) + 



n 



sinh T] 



:i-u' 



Vl>. = / dr] C4a27i^ 
Jo 



m^ cosh T] sinh 77 + n^ ( 1 + cosh 77) (coth 77 — 1 ^ 



.^2 



[u 



(30) 



t/33 = t/3-3 + Ul, 



1 I A 2 s^ 2 3 coth ?7 — 1 

dri< ApO ZTT —2 (± — u ; -r <-^4t 

sinh ?7 



:i - u') + C4a47r'rM 



2/„./n2 



(31) 



U,, = U^, + [/?,, [/f, = /;- rfr/ ^2 ^3^2 3 coth r/ 1 ^^ ^ ^,^^ 



f/* = / dry C4a27r^ 



sinh 7] 



1-^2 



+ (1 - u')Vm,n{v) 



(32) 



P^12 =W^f2 + W^2 

W^2 = - Apa^TT^ 



m. / (i?7< cosh(m?7) — sinh(|m|?7) >(coth?7 + \m\) 



X coth?7 



u 



Uy/l — M^' 



Vl — M^ sinh ?7 



f°° , [3 coth ?7 I 1 3 coth^ 77 — 2 2 coth 77 

+ m d7]< . . + \m\ —-2 + ^ . ,2 

Jo [ smh ?7 smh rj smh 77 



X < cosh(m?7) — sinh(|m|?7) > sinh t]uV1 — v? 



Wt = - C4a2n' 



- — - / (i?7< cosh(m?7) — sinh(|m|?7) 
m Jo I 



sinh 7] cosh 77^^ ^(r/)^' vl 



u^ 



{u'fu 



VT 



u^ 



+ m (i?7< cosh(?7i?7) — sinh(|?7i|?7) > sinh ?7 cosh; 



{u'fu 



(33) 



The conjugate momenta corresponding to Xj, A and B are Pj := ^^, /fj : = 
^^, Lj := ^^, since the variables a^, 6^ defined by (J2S1) are regarded as generalized 
velocities in the terminology of analytical dynamics. Here L and K are body-fixed spin 



and body-fixed isospin angular momenta. (Isospin K becomes the angular momentum 
operator acting on the target space.) The Hamiltonian defined through the Legendre 
transform, H^ '■= PjXj + Kjttj + Ljhj — Lp, is expressed in terms of Kj, Lj, Pj after 
eliminating Xj, aj, h^. Here it should be remarked that there is a primary constraint: 



which yields a selection rule for possible glueball spectrum below. 
The Hamiltonian in the rest frame Xj = reads 



(34) 



1 



Hf =M + - 



+ 



2t/n 
1 



w?. 



{Kl + Kl] 



1 



1 



Vii 



2^11 



-{Ll 



W^ 



12 



Wl2 



-K1L2 + 



W, 



21 



UuVn 
W21 






Ll) 



2U- 



-Ki 



(35) 



33 



Note that Uu = U22 diverges Un = U22 = OiV) ^ 00 as the volume V of the 
three-dimensional space. Therefore the Hamiltonian reduces to 



Hv 



M + 



21^11 



(Ll 



Ll) + 



2t/33^ 



M + 



2V, 



11 







(36) 
. Thus we obtain the 

(37) 



The Hamiltonian has the same form as the symmetrical top 
energy eigenvalue in the rest frame of the knot soliton: 



E 



where we have used the relationship L^ = J^ between the body-fixed spin angular 
momentum L and the space-fixed (coordinate-fixed) spin angular momentum J, since 
they are related to each other by spacial rotations: Jj = —Rjk{B)^Lk. For rotations 
in the target space, only rotations around the third axis are compatible with the 
boundary conditions and hence we have a relationship K^ = /| between the body- 
fixed isospin angular momentum K and the space-fixed (coordinate-fixed) isospin 
angular momentum / which are related to each other by rotations: Ij = —Rjk{A)Kk. 



Thus the energy of the knot soliton is modified after quantization into 



E = AfC^'M, + ApQ"'/' 



1 



n 



2V,\ 



'^(^ + l) + o TTT-TTrK: 



2 \U. 



33 



V, 



11. 



(38) 



where we have used the scaling argument for inertial tensors fj30|) -()32 |) and a must be 
replaced by a in V^^ and f/33, which follows from the substitution of 



4 Symmetries and constraints 

Every symmetry of a classical configuration induces a loop in configuration space. 
After quantization, these loops give rise to constraints on the wavefunction. In the 
foUowings, we study symmetries given by a rotation by a in space and by a rotation (3 



9 



in target space (isorotation). For example, one can consider a special path associated 
with the symmetries parameterized by a angle about k axis in space and /3 angle 
about j axis in target space: 

A(/5) := e'^^iAe-'^^^ = Ae''^"^!'^, B{a) := e'""^^ Be-'"^" = e-^""^/^^, (39) 

which follow from the commutation relations 

[K^, A] = ^Aa,, [Lfc, B] = -^a^B. (40) 

Then the end-point operator is given by 

L{a,p)=e"^^^e'''^K (41) 

If such rotations become symmetries inducing loops in configuration space, the Finkelstein- 
Ruhinstein (FR) constraint ^Hj niust be imposed on the wavefunction ip after quan- 
tization: 

— * — * — * 

exp(ipN ■ K) exp(2ae ■ L)ip = Xfr'^P, (42) 

where eTis the direction of the rotation axis in space, A^ is the rotation axis in target 
space (isospace), L and K are the angular momentum operators in space and target 
space, respectively. 

A scalar field theory can be quantized by considering the wave function on con- 
figuration space. This gives a purely bosonic theory. If the fundamental group of the 
configuration space is Z2, one can also consider wave functions on the covering space 
of configuration space. This gives rise to a quantum theory which contains fermions, 
provided so-called Finkelstein-Rubinstein constraints are introduced. 

In the quantization of solitons, therefore, two choices are possible; fermionic quan- 
tization or bosonic quantization. In the fermionic quantization, the FR phase xfr is 
given by 

{-|-1 if the induced closed loop is contractible 
— 1 otherwise (non-contractible) 

Whereas the bosonic quantization adopts the FR phase 

XFR = +1, (44) 

irrespective of whether the induced closed loop is contractible or not. 

There are several arguments why the Skyrme model has to be quantized fermion- 
ically. For a rotation by 27r, the following facts are known for the Skyrme model. 
Williams [201 and Giulini [2^ have shown that the configuration space admits spino- 
rial states, i.e., non-contractible loop with xfr = — 1 in the sectors of odd topological 
charge (winding number or baryon number) B, and no spinorial states, namely, only 
contractible loops with xfr = +1 are induced in the sectors of even topological charge 
B. Finkelstein and Rubinstein ^H] have shown that a 27r rotation of a Skyrmion of 
charge B is homotopic to an exchange of two Skyrmions of charge B. This implies 
that an exchange of two identical Skyrmions gives rise to xfr = — 1 if and only if 
their topological charge B is odd, or xfr. = +1 if and only if B is even. The Skyrme 
model allows for fermionic quantization with half-odd angular momentum in the odd 
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B sectors, and bosonic quantization with integer angular momentum only in the even 
B sectors. This is the topological spin-statistics theorem for the Skyrmion [22] • Kr- 
usch j2n] has shown that a 27r isorotation of a Skyrmion also gives rise to Xfr = — 1 
if and only if B is odd. 

For Skrymions which are well-approximated by the rational map ansatz (rational 
map Skyrmion), a simple formula for the FR phase in the fermionic quantization has 
been given j2ni, Xfr = (—1)^, N = B{Ba — /3)/27r, and the formula was generalized 
to more complicated situations in ^^. 

Here it is instructive to recall the consideration done by Braaten and Carson |T7j 
for the B = 2 Skyrmion under the product Ansatz. For the B = 2 Skyrmion, 
the paths induced by vr rotations about coordinate axes are closed in the sense of 
returning the transformed field at 6 = n to its original value at ^ = 0: U{f, t) = 
exp{ii3N ■ K) exp{iae ■ L)U{f,t)exp{—iae- L) exp{—ij3N ■ K) for the Skyrme field 
U{f,t). Consequently, acting the end-point operator exp{ipN ■ K) exp{iae ■ L) on 
the wavefunction gives rise to a phase xfr = +1 or —1 depending on whether the 
associated path is contractible or not, respectively. The analysis of fTj shows that 
Tc rotations about x and y axes with end-point operators e^^^i^^^^i and e^'^^ie^'^^z 
correspond to non-contractible paths leading to Xfr = ~1; while vr rotation about 
the z axis with the end-point operator e*'^^^ corresponds to a contractible path leading 
to Xfr = +1- 

In this paper, we intend to use the Faddeev model for describing effectively glue- 
balls which are bosons with integer spins to be identified with composite particles 
made of spin-one gluons. In this paper, therefore, we choose the bosonic (collective 
coordinate) quantization of knot solitons in the Faddeev model. (This choice should 
be compared with the fermionic quantization chosen in Krusch and Speight JH]-) In 
view of this, we study the symmetries and the corresponding constraints on the wave- 
function to be imposed after quantization. In the Faddeev model under the simple 
toroidal ansatz, the minus sign in (fTTj) or (fTH|) can not be expressed by using A{i3) and 
B{a) in the manner described in ()39|1 . Rather, it denotes a reflection in the target 
space. In any case, therefore, the spacial rotations by ir radians about coordinate axes 
X and y do not represent closed loops in the configuration space of the color field, 
while the spacial rotation by vr radians about z axis (fTIHl represents a closed loop, as 
a special case of the continuous axial symmetry if —^ (f + 'j, i.e., 7 = vr. Therefore, to 
obtain closed loops induced by rotations about coordinate axes x and y, the rotation 
angle must be at least 2vr radians. 

In quantum system, the operators J^ = L'^jL^, J3, I^ and K^ form a set of com- 
muting observables and the eigenfunction ip can be labelled by the quantum numbers 
L,L3, Js,!^ and K3. We consider states with given J and I3, namely, L = J and 
K3 = iJa and L3 = —nK^. Then the eigenfunction for the above Hamiltonian is given 
up to a normalization constant hy ip = | J, Is, K^) x |L, J3, L3) where | J, /s, K^) corre- 
sponds to a finite rotation in the target space, i.e, isorotation A G SU{2)/Z2 = 5*0(3), 
and I J, J3, L3) to a finite spacial rotation B G 5*0(3). 

Therefore, reflecting the transformation law ()17p. ()18|) and p9p of the classical 
configuration for the vr rotations about x, y and z axes, the end-point operators 
representing closed loops, namely, the 2vr rotations about x and y axes and the vr 
rotation about z axis are constructed to act on the wavefunction tp = \I,I^,K-i) x 
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|L, J3, L3) in the following way: for n = even (L3 = —nK^ = integer): 

^2.^K,^2.^L,^ =(-1)2^+2^^^ = ^^^^, (45a) 

While, for n = odd (L3 + Ji'3 = (1 — n)i^3 = integer), we have 

^2.iK,^2.iL.^ ^(_l)2L+2/^ ^ ^^^^^ (4g^) 

^2niK^^2niL,^ =(-1)2^+27^ ^ ^^^^^ ^4g^) 

^^iKs^^iU^ =(^_l^L,+K,^ ^ (_l)(l-")i^3^ = ^pj^^_ (46c) 

Here we have used the properties of the Wigner D-function PH]: {J, M'|e*'^"^^| J, M) = 

SM',~M{-iy^^^ and (J,M'|e^"-^=^|J,M) = 5A//,_M(-l)'^+*'^e-'"^, leading to e*"-^- 1 J, M) 
(-l)-^+^^e-*^^V, -^) and e*^'^^|J,M) = {-iy+^^J,-M), in addition to the rela- 
tion e^'^-^^l J, M) = e^'^^l J, M) = (-1)*V> M). They are applied to our case to obtain 
the formulae: 

e'^^^m J3, Ls) =(-l)^+^^e-^'^^«|L, J,, -L3), 

e^'^^^lL, J3, L3) =(-l)^+^=^|L, J3, -i^s), 

e^'^^^lL, J3,L3) =(-l)^«|L, J3,L3), (47) 

and 

Besides rotations considered in the above, closed loops are indeed obtained by 27i 
rotations in space and the target space, and our restriction for the FR constraint for 
these loops in the bosonicc quantization xfr = +1 leads to 

e2-^3^ = (-1)'^^^ = XfrV', e^-^^^^ = (-1)'^^^ = XfrV', (49) 

which implies that L3 and K^ must be integers. 

Moreover, we examine the restriction on quantum numbers coming from the FR 
constraints for the symmetries (jl^j) and ^^ in the bosonic quantization xfr = +1- 
A condition, I + L = integer, must be satisfied for any n. For n = even, L3 = —nK^= 
even integer and hence L = J must be integer (excluding half-odd integer) yielding 
that I is also an integer. For n = odd, we have also a condition, I+L = integer. In this 
case, we can conclude only that L3 and K3 are integers. Hence L = J is also an integer 
yielding that / is also an integer. Thus, both I and J are integers. These results are 
reasonable for our purposes of identifying the Hopf solitons with bosons and can be 
a self-consistency check for choosing the bosonic quantization. Further symmetries 
may impose additional constraints for quantum numbers, leading to further selection 
rules for possible spectrum of glueballs. 
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SU(2) Lattice gauge theory 


J^^' 


Mg/^ 


Mg/Mo++ 


0++ 


3.74 ± 0.12 


1 


2++ 


5.62 ± 0.26 


1.46 ± 0.094 


0-+ 


6.53 ± 0.56 


1.78 ± 0.24 


2-+ 


7.46 ± 0.50 


2.03 ± 0.20 


1++ 


10.2 ± 0.5 


2.75 ± 0.15 


1-+ 


[10.4 ± 0.7] 


[3.03 ± 0.31] 


3++ 


9.0 ± 0.7 


2.46 ± 0.23 


3-+ 


[9.8 ± 1.4] 


[2.91 ± 0.47] 



Table 1: Glueball masses, (right column) in units of the lightest scalar glueball mass 
and (left column) of the string tension, in the continuum limit of SU(2) Yang-Mills 
theory on a lattice J2S|- Values in brackets have been obtained by extrapolating from 
only two lattice values and so should be treated with caution. 



5 Glueball mass, gluon mass and vacuum condensate 

As in the previous section, we restrict the following considerations to bosonic quanti- 
zation with the trivial FR phase Xfr = +1- Then the eigenfunction has the simplified 
form: 



ip = \L,Ls,K3), 

J = 0,l,2,..., \Ks\ 



0,1, 



J, 



0,±1, 



±J, 



(50) 



where we have suppressed J3 and Is, since we find no constraints on the values of 
J3 and /s according to ^Hl- This wavefunction is identified with the Wigner D- 
function ^H] i^LsK-i ~ ^^LiKai^^P^l) f*^^ ^ finite rotation with the Euler angles 



a, /5, 7 where C is a normalization constant C = J {2 J + l)/(87r2), which is known as 
the eigenfunction of the above Hamiltonian (jHBj) of a symmetrical top. 

In the followings, we consider only quantum states with Hopf charge \Qh\ = 
1,2 as candidates for glueball states constructed from gluons with spin-one. This is 
because the simplified toroidal ansatz adopted in this paper is valid only for \Qh\ = 
1,2 according to numerical calculations ^21- See Tabled for results of numerical 
simulations for glueballs on a lattice. For \Qh\ = 1, we have a case {n,m) = (1, 1). 
For \Qh\ = 2, our numerical calculations for Qh = 2 show that the classical energy 
M* (the rest mass of the knot soliton) for {n,m) = (2, 1) is lower than the value for 
{n,m) = (1,2), see Table El 

For m = n, i.e., Qh = m"^ = 1, 4, 9, ■ ■ ■ , the Vakulenko-Kapitanskii lower-bound 
on the energy [211 by the Hopf index, E > const. [Qj:/!^/^, is saturated ^7\ by the 
classical background configuration, although the generic Hopf soliton is not of the 
BPS type. In this case, the lowest 0"*"+ glueball mass is expressed as 



M 



o+^ 



1 /2 

const. AfC4 \Qh 



1 3/4 



const. (? 



-1 



-K)\Q 



',H 



1 3/4 



(51) 



by the energy of the Hopfion with non-zero Hopf index Qh- The non-perturbative 
dependence g~^ of the mass on the coupling constant is characteristic of the soliton. 
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SU(2) Faddeev model 


Qh 


{n,m) 


M, 


V,\ 


uh 


W^l*2 


a 


1 


(1,1) 


292.1 


1073 


1237 


148.3 


1.177 


2 


(2,1) 


465.3 


1274 


432.8 





1.740 


2 


(1,2) 


572.6 


3049 


2550 


941.2 


1.310 


3 


(3,1) 


647.7 


2407 


416.6 





2.293 


3 


(1,3) 


917.2 


7357 


5246 


900.5 


1.416 


4 


(2,2) 


845.3 


3260 


977.5 





1.756 


4 


(4,1) 


860.3 


5317 


431.1 





3.128 


4 


(1,4) 


1312 


20454 


17974 


2007 


1.498 



Table 2: The (dimensionless) mass M* and the size a of the knot sohton and non- 
vanishing and finite components V^\, U^^, W12 of the inertial tensor for various Hopf 
indices Qh = rnn at Ap = 1 and C4 = 1 for the Faddeev modeL 



For definiteness, we restrict the following consideration to the case (n, m) = (2, 1), 
since the other cases can be treated exactly in the same way. We adopt a simple 
identification as inputs: 



0++ ^ \L, L3, Ks) = |0, 0, 0) (m = 1, n = 2; Qh = 2), 
2++ ^ \L, L3, Ks) = |2, 0, 0) (m = 1, n = 2; Qh = 2). 



(52) 



Note that Parity is given by P = ( — 1)"' for K^ = and the wavefunction reads 
Dqq'{(P,9, —if) = (^00(6') = Pj{cos9) with no yj-dependence for K^ = = L^ . Under 
this identification, the lowest two glueball masses are given by 



Mo+H 



AfCI^'m,, M2++ = AfCI^^M, + SApCf^^ 



/V, 



11, 



(53) 



which yield the mass ratio M2++/M0++ = 1 + 30^:'^ / {]VhV{^) . 

Our numerical calculations show M^ = 465, V{i = 1274, U^^ = 433 for [m = 
l,n = 2]Qh = 2). By combining these results with those of numerical simulations 
of Yang- Mills theory on a lattice [53]: M2++/M0++ = 1.46 ± 0.09, a parameter is 
determined to be C4 = [3/{M^V{^{M2++/Mo++ - l))]!/^. 

C4 = 0.00303 ~ 0.00369. (54) 

The ratio of the lowest glueball mass to the gauge-invariant gluon mass [2j due to ^i 

V(_x2)ym = AfC4~'/', 



M 



X 



is determined from the ratio M2++/M0++ alone: 

Mo++/Mx = C4M, = 1.41 ~ 1.72. 



(55) 



(56) 



This results supports that a glueball is composed of gluons. Another parameter Ap 
is determined from an input mass Mq++ = (3.74 ± 0.12)^/0- ~ 1.5GeV to be 



Ap = Mo++/{Cl^^M,) = 0.0530 ~ 0.0586GeV. 



(57) 
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SU(2) Faddeev model 


Qh 


{n,m) 


C4 


QM, = Mo++/Mx 


AF[GeV] 


a[fm] 


Mx[GeY] 


1 


(1,1) 


0.00456 


1.332 


0.0760 


0.206 


1.126 


2 


(2,1) 


0.00332 


1.543 


0.0560 


0.352 


0.972 


2 


(1,2) 


0.00193 


1.107 


0.0596 


0.190 


1.356 



Table 3: The parameters of the Faddeev model determined from the identification 
with the glueball masses. Only the central values are indicated. 



The radius of the torus is a = Ap^C^^a = 1.64 ~ 1.99 (GeV)"^ = 0.32 ~ 0.39 fm, 
which is the associated Compton wavelength to be compared to that of the pion 
with mass m^ = 140 MeV, 0.197/0.14 = 1.41 fm. This result is consistent with the 
prediction of phenomenological approaches, see e.g., j^Hl- Thus we obtain an estimate 
of the gluon mass 



Mx = 0.87 ~ 1.41 GeV, 
and the numerical value of the dimension two condensate: 



i-K) 



YM 



A2 = (0.0530 ~ 0.058GeV)2. 



(5J 



(59) 



All the cases considered are summarized in Table El The resulting values for Mx 
are consistent with the result of numerical simulations on a lattice [SU] : 



M 



X 



1.2 GeV. 



(60) 



6 Conclusion and discussion 

We have performed the collective coordinate quantization of the classical knot soliton 
solution obtained under the toroidal ansatz of the Faddeev model. Then we have 
identified the resulting energy spectrum with the glueball mass spectrum of the SU(2) 
Yang-Mills theory, supposing the Faddeev model is a low-energy effective theory of 
SU(2) Yang-Mills theory. We have used two input masses Mo++ and M2++ obtained by 
numerical simulations of Yang-Mills theory on a lattice to determine two parameters 
in the Faddeev model for predicting glueball masses other than O"*""*" and 2"*"+. 

Thanks to the NLCV, we can expect that the gauge-invariant mass term ^M^^X^ 
can be induced in Yang-Mills theory, since the dimension two composite operator X^ 
can be defined in the gauge-invariant way based on the NLCV. From the collective 
coordinate quantization of the Faddeev model and the identification of the glueball 
mass spectrum, we have obtained the reasonable value for the gauge-invariant gluon 

X?, ) . These results are consistent with an observation that the 

A'/YM 

X2\ 



mass Mx oc \l g 

gauge- invariant vacuum condensation (^t) of mass-dimension two takes place in 

Yang-Mills theory and gives a gluon mass Mx oc 



9' 



X^, ) . The resulting mass 



scale enables us to derive the kinetic term in the Faddeev model as a LEFT of Yang- 
Mills theory, which is needed to avoid the no-go theorem of Derrick and to allow the 
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existence of the topological soliton as a result of balancing the kinetic term with the 
Skyrme term of four derivatives. Thus the dimension two condensate (X^) in QGD 
leads to the Faddeev model and the Faddeev model could be a realistic low-energy 
effective theory of SU(2) Yang-Mills theory. 

It has been also shown [2H1 that magnetic monopoles in QGD can be defined in 
the gauge-invariant way based on the NLCV. Then the Wilson loop average can be 
calculated from the linking number between the magnetic glue-knot and the electric 
Wilson loop. This will shed more light on the existence of mass gap in QGD and 
quark confinement in QGD. Extending this work to color SU(3) case will be done in 
a subsequent paper. 
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